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Looking for the Top-squark at the Tevatron with four jets 
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The scalar partner of the top quark is relatively light in many models of supersymmetry breaking. 
We study the production of top squarks (stops) at the Tevatron collider and their subsequent decay 
through baryon-number violating couplings such that the final state contains no leptons. Performing 
a detector-level analysis, we demonstrate that, even in the absence of leptons or missing energy, stop 
masses upto 210 GeV/c^ can be accessible at the Tevatron. 



o 
o 

(N 



> 

ON . 

o ■ 
O ' 

00 

o ; 

in • 

o : 
i!^ ■ 

Oh! 
I 

Oh 



X 



PACS numbers: 14.80.Ly, 12.60.Jv 

The conservation of either of baryon (B) and lepton 
(L) numbers is not dictated by any fundamental principle 
and is but an accidental feature of the perturbative sector 
of the Standard Model (SM). Indeed, any explanation of 
the observed baryon asymmetry of the universe needs at 
least one of i? or L to be broken to a significant degree. 
Many extensions of the SM, whether supersymmetric or 
not, naturally admit both B and L violation and care 
must be taken that both are not violated strongly so as to 
render the proton very unstable. As can be easily appre- 
ciated, breaking of either B or L would lead to significant 
alteration in phenomenology, and in particular, collider 
signatures for physics beyond the SM. While numerous 
studies have been undertaken in the context of L viola- 
tion, in this article we seek to examine the experimentally 
more challenging case of a broken baryon-number. 

Low energy supersymmetry (susy) is widely consid- 
ered to be a benchmark for scenarios going beyond the 
SM. Since the most general renormalizable Lagrangian 
consistent with both gauge symmetries and SUSY contain 
terms that break both B and L, stability of the proton 
is normally ensured by the imposition of an ad hoc dis- 
crete symmetry, namely i?-parity 1]. However, since the 
same end can be achieved by the imposition of L alone, 
we allow, in the superpotential, terms of the form 



(1) 



where C/]j and D^j^ denote the right-handed up-quark and 
down-quark superfields respectively. The Yukawa cou- 
plings A"^j, are antisymmetric under the exchange of the 
last two indices. The corresponding Lagrangian can then 
be written in terms of the component fields as: 



A 



u'idUl+uld^Al]+\i.c., (2) 



thus allowing a squark to decay into a pair of quarks. 
While resonanct production in a hadron collider is pos- 
sible as well 0, the corresponding rates can be appre- 
ciable only if two of the superfields belong to the first 
generation, and then too are limited by the size of the 
couplings A". As can be expected, the latter are con- 
strained by various low-energy observables HH, though 



the couplings involving the second and third-generation 
fields alone can be relatively large |S] • It is thus advisable 
to concentrate on the (model-independent) strong inter- 
actions for squark production and consider the effect of 
A" only in the decays. 

In most SUSY models, the large top Yukawa coupling 
results in the the Hghter stop, ti, being hght compared 
to the other squarks. Since the realization of the mech- 
anism of electroweak baryogenesis within the context of 
the MSSM requires light stops, with masses of about or 
smaller than the top quark mass 0, there is an added 
motivation to consider such scenarios. 

At hadron colliders, stop production proceeds over- 
whelmingly via the strong interaction and the corre- 
sponding cross sections are well known at leading order 
ij. The next-to-leading order QCD and SUSY-QCD cor- 
rections have been computed |3 and implemented nu- 
merically in PROSPINO [y, 01 , which we use along with 
the CTEQ5 parton distribution functions JTj. We fur- 
ther assume that the masses of the gluino and the other 
squarks are larger than about 250 GeV so that they do 
not alter the NLO cross section significantly 0. This, fur- 
thermore, precludes any significant enhancement of the 
stop production cross sections via cascade decays thereby 
making our estimates conservative. 

The prospects for stop discovery at the Tevatron have 
been examined both in the context of i?-conserving su- 
pergravity inspired scenarios |12| as well as in the con- 
text of low-energy SUSY breaking^^. Search efforts at 
the LEP and the Tevatron, irrespective of the stop decay 
mode, have only proved unsuccessful Q. The reach, 
at Run II, depends crucially on the decay chain (and, 
hence, the SUSY spectrum) and, for an integrated lumi- 
nosity of 2 fb~^ typically ranges between 165-190 GeV. 
For a stop light enough such that ti — > x+6 is kinemat- 
ically forbidden, the details of the decay depend very 
sensitively on the mass splitting between the stop and 
the lightest neutralino (note that if i?-parity is broken, 
the stop is even allowed to be the lightest SUSY parti- 
cle). If, for example, < mw + to& + w^Oj only two 
_R-conserving decay modes are kinematically accessible, 
namely (i) the loop- induced flavor-changing two-body de- 
cay t cx\ and (m) the four-body decay via a virtual 
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W boson, i or ivbxi- It is easy to 

see that either of these partial widths are small and may 
be superseded by i?-violating modes even for moderate 
values of A". For the rest of this paper, we shall assume 
that at least one of the modes ti ^ b + s (d) has a sig- 
nificant branching fraction. (We refrain from discussing 
ti — > d + s for reasons of experimental sensitivity.) 

At this stage, we digress to point out that the stop (or 
any other squark) is not the only conjectured strongly- 
interacting particle that may decay into a pair of quarks. 
Even in the simplest nonsupersymmetric grand unified 
theories (GUTs), B may be violated in both the gauge 
and the scalar sector interactions. The corresponding el- 
ementary particles, namely diquarks, can be either spin-0 
or spin-1 and have baryon and lepton numbers 2/3 and 
respectively Tsj. A generic diquark may transform as 
3 or 6 under SU{3)c, as triplet or singlet under SU{2)l 
and can have electric charges \Qv\ = 1/3,2/3 or 4/3. 
Compared to the A"s, diquark couplings are typically 
less restricted both in terms of symmetry requirements 
(allowing, for example, the experimentally easier mode 
V ^ b + b 0]) as well as low-energy constraints As 
far as scalar diquarks are concerned, a SU (3)c triplet has 
the same production cross section (and phase space dis- 
tributions) as a stop of identical mass, while a sextet has 
a larger one on account of the color-factor. The cross 
section for a vector diquark depends on the exact nature 
of its gauge interactions and is significantly larger. More- 
over, a generic diquark tends to decay dominantly into a 
pair of quarks. The stop, thus, is the most conservative 
choice from this genre. 

At the partonic level, our final state, thus, consists of 
(b q) (bq) where q is either a d- or a s-quark and the paren- 
thetical pairing is to denote that the combinations arise 
from the decay of an (anti-)stop. The SM backgrounds 
were generated with both madgraph 17] and PYTHIA 
6.206(131 and tested for consistency. Using the latter, we 
generate complete events with initial and final state radi- 
ation, multiple interactions, etc., and complete evolution 
(hadronisation and decays) of the partons into final state 
particles. The latter are passed through a toy detector 
simulation (using tools in PYTHIA) and event recon- 
struction algorithm mimicking a typical Tevatron Runll 
detector. The toy calorimeter has cell sizes of Arj = 0.1 
and A(j) =15°. Jet reconstruction has been done employ- 
ing the cone algorithm with AR = yj (Ar;)^ -I- (Ac/))^ < 
0.7 and using calorimeter clusters with Et > 1-0 GeV as 
seeds for jet formation. Only jets with \r]jet\ < 2.4 and 
Et > 15 GeV and leptons with \rii\ < 5 are considered. 
Tagging of fe-jets has been done using decay lengths of 
&-hadrons such that ~ 60% of tt events have at least one 
6-jet tagged[li|. Apart from vertex tagging, soft lepton 
tag may be used to enhance 6-tagging. Since the event 
features used in this analysis, viz., jet and lepton Px and 
ry, jet multiplicity, I^t and &-tag are rather robust and 
easy to implement, our results would be fairly indepen- 



dent of the detailed features of a particular detector. 
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FIG. 1: Differential distributions (normalized to unity) in 
different variables for both signal (m^^ = 120 GeV/c^) and 
various SM backgrounds. The top two panels correspond only 
to the acceptance cuts. For the bottom one, all selection 
critera of eq.J^I other than that on TVjct have been imposed 
as well. 

The signal events would be characterised by four hard 
jets, two of them being 5-jets. Since leptons or neutri- 
nos in such events can occur only as decay products of 
hadrons, these would be soft. This then inspires our se- 
lection criteria 

< 15.0 GeV/c no lepton with > 15 GeV 



Ej^t Et > 200 GeV 



btag 



Mbb,Mjj ^ (70, 100) GeV/c^ 



(3) 

Apart from the ZZ process (which is largely elimi- 
nated by the last requirement above), backgrounds also 
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arise from ti events with both tops decaying hadronically 
(these typically have more than four jets, see Fig. ^ as 
well as bb events accompanied by either or both of mul- 
tiple interactions and hard gluon or photon radiation. 
Although bb events have a huge cross-section, the cut on 
J2j Eip is very effective with a rejection factor of about 

10^ as Fig. ^amply demonstrates. 

This still leaves a large background. However, in 
tit\ events, of the two jet pairings viz. (61 ji, 62^2) 
and (&1J2, ^2^1 ), the one representing the decaying stops 
should be associated with only a small difference in the 
reconstructed invariant masses. Hence, our final selection 
criterion is that 



A^^ l < 20 GeV/c^ 



(4) 



for at least one pairing. For the signal events, the corre- 
sponding average of the two masses is expected to show 
a sharp peak around m^^ as is evinced by Fig. |21 whereas 
the other pairing has a rather flat distribution. 
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FIG. 2: Distribution of the average reconstructed mass for 
tit* events for m^^ = 120 GeV/c? (10^ events generated). 
The dark (black) line corresponds to the combination with 
the smaller difference between the two invariant masses; the 
light (purple) line represents the other combination. 

We have simulated 10^ events for each mr and also for 

_ II 

tt and ZZ events. Though bb events have a very small 
selection efficiency, they have a very large cross-section, 
and constitute the bulk of the background events passing 
the selection cuts. Hence, a very large set of bb events 
2.5 X 10*) have been generated to get a good estimate of 
the background distribution. As for the signal events, for 
low m^^ , a large fraction of the events fail to satisfy the jet 
selection criteria leading to a small selection efficiency e 
(FiglSJ. As m^_^ increases, the situation improves rapidly; 
however beyond 150 GeV/c^, this effect saturates and is 
more than offset by the rejection on account of hardening 
of lepton Pt and I^t ■ The rapid fall in the effective cross- 



section (cr ■ e) is, of course, reflective of the p-wave nature 
of scalar production. 
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FIG. 3: (a) The detection eficiency for stop-pair production 
and its product with the cross-section as a function of the 
mass, (b) The distribution in the average of the two masses 
Mhj corresponding to the minimum difference. While the 
signal profile for three different stop masses are given by the 
points, the solid and dashed lines show the backgrounds from 
bb (overwhelmingly dominant) and tt events. The ZZ rate 
falls below the scale of the figure. 

Whereas the signal events show a sharp peak in M^J^, 
the background is much flatter (Fig. This allows us to 
identify a range in M^J^ where the signal is most signif- 
icant and calculate the x^- Working with a conservative 
choice of a 50 GeV/c^ bin, we use this to obtain an 
exclusion plot in the BR{ti — > bq) — mj_^ plane (Fig 0)) 
that may be reached by the Tevatron experiments. With 
as little as 2 fb~^ data, such an analysis would have a 
reach upto 185 GeV/c^ (for BR(ti bq) = 100%), and 
on the other hand probe down to BR{ti — > bq) ~ 4% for 
TOj = 70 GeV I (? . Similarly, we may be able to put an 



• bq) for stop masses upto 
^. A combined analysis of 



upper bound on the BR(ti - 
200 GeV/c^ with dnt = 4 fb 
the data from the two Tevatron Runll experiments would 
serve to push the limits even further further. 

In summary, we have outlined above a robust stop- 
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search strategy based on selection criteria which are easy 
to implement. Most importantly, a final state devoid of 
leptons or missing energy is shown to be very promis- 
ing and competitive with other modes that have been 
used so far. Specific features dependent on detector capa- 
bilities may be used, particularly in a multidimensional 
analysis, to better discriminate signal and background 
and probe added regions in the parameter space. Fur- 
thermore, the sensitivity can be enhanced by considering 
TOjj^ -dependent selection criteria rather than the univer- 
sal cuts that we have chosen to impose. In fact, even an 
analysis of the currently accumulated data would serve 
to probe a significant region of the SUSY parameter space 
that has not lent itself to an examination so far. And as 
we have already pointed out, the analysis is not limited to 
the top-squark or supersymmetry alone but can be read- 
ily extended to diquarks, which, in fact, are generically 
associated with even larger cross sections. 



FIG. 4: Exclusion contours at 90% CL in the BR{ti bq) - 
m^^ plane that may be achieved for different values of total 
integrated luminosity. 



DC would like to thank the DST, India for financial as- 
sistance under the Swarnajayanti Fellowship grant. MM 
would like to thank S. Chakraborty for useful discussions 
and TIFR for the use of computing facilities. 



[1] P. Fayet, Phys. Lett. B69 (1977) 489; G. Farrar and 
P. Fayet, Phys. Lett. B76 (1978) 575. 

[2] Such scenarios can be motivated from a class of super- 
symmetric GUTs as well, e.g. K. Tamvakis, Phys. Lett. 
B382 (1996) 251. 

[3] F. Abe et al. (CDF CoUab.), Phys. Rev. D 55, 5263 
(1997). 

[4] J.L. Goity and M. Sher, Phys. Lett. B346 (1995) 69; 
C.E. Carlson, P. Roy and M. Sher, Phys. Lett. B357 
(1995) 99; B. Brahmachari and P. Roy, Phys. Rev. D 50, 
R39 (1995). 

[5] G. Bhattacharyya, D. Choudhury and K. Sridhar, Phys. 
Lett. B355 (1995) 193; D. Chakraverty and D. Choud- 
hury, Phys. Rev. D 63, 075009 (2001); ihid.Qi, 112002 
(2001); 

[6] G. Bhattacharyya, arXiv:hepph/9709395; R. Barbier et 
al, arXiv:hep-ph/9810232 B. C. Allanach, A. Dedes and 
H. K. Dreiner, Phys. Rev. D 60, 075014 (1999); R. Bar- 
bier et a/., arXiv:hep-ph/0406039 

[7] See, for example, C. Balazs, M. Carena and 
C. E. M. Wagner, Phys. Rev. D 70, 015007 (2004) and 
references therein. 

[8] P. R. Harrison and C. H. Llewellyn Smith, Nucl. Phys. B 
213, 223 (1983) [Erratum B 223, 542 (1983)]; S. Dawson, 
E. Eichten and C. Quigg, Phys. Rev. D 31, 1581 (1985). 

[9] W. Beenakker et al, Nucl. Phys. B 515, 3 (1998). 



[10] W. Beenakker, R. Hopker and M. Spira, 

arXiv:hep-ph/961 1232 
[11] H. L. Lai et al [CTEQ Collab.], Eur. Phys. J. C 12, 375 

(2000). 

[12] H. Baer et al, Phys. Rev. D 44, 725 (1991); H. Baer, 
J. Sender and X. Tata, Phys. Rev. D 50, 4517 
(1994); R. Demina, J. D. Lykken, K. T. Matchev 
and A. Nomerotski, Phys. Rev. D 62, 035011 (2000); 
A. Djouadi, M. Guchait and Y. Mambrini, Phys. Rev. D 
64, 095014 (2001). 

[13] S. Ambrosanio et al Phys. Rev. D 54, 5395 (1996); 
C. L. Chou and M. E. Peskin, Phys. Rev. D 61, 055004 
(2000); R. Culbertson et al, arXiv:hep-ph/0008070| 
M. Carena et al, Phys. Rev. D 66, 115010 (2002). 

[14] See, S.P. Das et al, Phys. Lett. B 596 293 (2004) and 
references therein. 

[15] F. Zwirner, Int. J. Mod. Phys. A3 (1988) 49; J.L. Hewett 
and T.G. Rizzo, Phys. Rep. 183 (1989) 193 and refer- 
ences therein. 

[16] This is currently under study. 

[17] F. Maltoni and T. Stelzer, JHEP 0302, 027 (2003). 
[18] T. Sjostrand et al, Comp. Phys. Comm. 135, (2001) 238. 
[19] D. Acosta et al (CDF Collab.), Phys. Rev., D71, 072005 

(2005); ihid. D71, 052003 (2004); V.M. Abazov et al (DO 

Collab . ) , arXiv:hep-ex/0504058 



